The stiffness exponent of two-dimensional Ising spin glasses for non-periodic 
boundary conditions using aspect-ratio scaling 
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We study the scaling behavior of domain-wall energies in two-dimensional Ising spin glasses with 
Gaussian and bimodal distributions of the interactions and different types of boundary conditions. 
The domain walls are generated by changing the boundary conditions at T = in one direction. 
The ground states of the original and perturbed system are calculated numerically by applying an 
efficient matching algorithm. Systems of size L x M with different aspect-ratios 1/8 < L/M < 
64 are considered. For Gaussian interactions, using the aspect-ratio scaling approach, we find a 
stiffness exponent 6 — —0.287(4), which is independent of the boundary conditions in contrast to 
earlier results. Furthermore, we find a scaling behavior of the domain-wall energy as predicted 
by the aspect-ratio approach. Finally, we show that this approach does not work for the bimodal 
distribution of interactions. 
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I. INTRODUCTION 



Although spin glassesEI have been studied for more than 
two decades, finite-dimensional spin glasses are still far 
from being well understood. This is true alaa-ituc? = 2 
dimensions, where no stable spin-glass phasenoDO exists 
at T > 0. The fact that the transition temperature, T c , 
is zero can be seen by studying domain walls between 
ground states induced by changing the boundary condi- 
tions. For a system of size L x M, if the domain wall is 
forced to run in the ^-direction, the mean domain-wall 
energy scales like M e in the thermodynamic limit, where 
9 is the stiffness exponent. IfO < 0, as in two-dimensional 
spin glasses, no stable phase can exist at finite tempera- 
ture. 

Recently!! slightly different values of the stiffness ex- 
ponent have been found for Gaussian distribution of the 
bonds for different choices of the boundary -conditions. 
Here, by applying the aspect-ratio methodB, we show 
that this difference is due to corrections to scaling which 
persisted, at least for one set of boundary conditions, 
even for the quite large system sizes used in Ref. ^. The 
results presented here are consistent with = —0.287(4) 
for all boundary conditions. This agrees with the result 
found in Ref. |^ where periodic boundary conditions were 
applied in all directions for smaller systems than investi- 
gated here, indicating a high degree of universality for the 
exponent 6. Apart from using larger systems and differ- 
ent boundary conditions, we go in another respect beyond 
the previous work, because aspect ratios, R = L/M, less 
than unity are studied here. This enables us to verify the 
proposed scaling behavior in the region R < 1, where the 
finite-size scaling depends on the boundary conditions. 

The model we study consists of N = L x M Ising spins 



Si = ±1 on a square lattice with the Hamiltonian 



i j Si Sj , 



(1) 



where the sum runs over all pairs of nearest neighbors 
and the Jij are quenched random variables. Here, 
we consider two kinds of disorder distributions: 

(i) Gaussian with zero mean and variance unity, 

(ii) bimodal, i.e. Jij = ±1 with equal probability. 

To study whether an ordered phase is stable at fi- 
nite teismefatures, the following procedure is usually 
appliedaclalZl. First a ground state of the system is calcu- 
lated, having energy Eq. Then the system is perturbed 
to introduce a domain wall and the new ground state en- 
ergy, i?Q Crtb is evaluated. For each sample, the domain- 
wall energy is given by -Edw = |£Q Crtb — Eq\. Here, we 
apply free boundary conditions (be) in the y-direction. 
To study the scaling behavior of the domain- wall energy, 
we consider two means to introduce domain walls: 

(i) "P-AP": First a ground state with periodic be in 
the x-direction is calculated. Then the system is 
perturbed by introducing antiperiodic be in that 
directions, e.g. by reversing one line of bonds par- 
allel to the y-direction. 

(ii) "F-DW" (called "F-AF" in Ref. |): First a ground 
state with free be in the x-direction is calculated. 
For the perturbed system, we then add extra bonds 
which wrap around the system in the x-direction, 
and have a sign and strength such that they force 
the spins they connect to have the opposite relative 
orientation to that which they had in the original 
ground state. 
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In both cases, M denotes the size of the edge along 
which the boundary conditions are changed to induce a 
domain wall (i.e. the size in the y-direction) . The scaling 
behavior 

(E DW ) = M e F{L/M) (2) 
has been predicted!, with the following limiting forms: 

{R- 1 for R -> oo 

R e-{d-i)/2 (p_ AP ) f or i?^o (3) 
R 6 -^- 1 ) (F-DW) for R -»• 0. 

We shall verfiy that the data does follow this scaling form 
very well. 



II. THE ALGORITHM 

In greater than two dimensions, or in the presence of a 
magnetic field, the exact calculation of spin-glass,|ground 
states belongs to the class of NP-hard problemsErcl. This 
means that only algorithms with exponentially increasing 
running times are known. However, for the special case 
of a planar system without magnetic field, e.g. a square 
lattice with periodic boundary conditions in at most one 
direction, there are efficient polynomial-time "matching" 
algorithmsEi. The basic idea is to represent each-realiza- 
tion of the disorder by its frustrated plaquettesO. Pairs 
of frustrated plaquettes are connected by paths in the 
lattice and the weight of a path is defined by the sum of 
the absolute values of the coupling constants which are 
crossed by the path. A ground state corresponds to the 
set of paths with minimum total weight, such that each 
frustrated plaquette is connected to exactly one other 
frustrated plaquette. This is called a minimum-weight 
perfect matching. The bonds which are crossed by paths 
connecting the frustrated plaquettes are unsatisfied in 
the ground state, and all other bonds are satisfied. 

For the calculation of the minimum-weight per- 
fect matching, efficient polynomial-time algorithms are 
availableqap Recently, an implementation has been 
presentecO, where ground-state energies of large systems 
of size N < 1800 2 wera-calculated. Here, an algorithm 
from the LEDA libraryuJ has been applied, which limits 
the system sizes due to the restricted size of the main 
memory of the computers which we used. Furthermore, 
each system border with a free boundary decreases the 
running time relative to a boundary with pbc. The rea- 
son is that frustrated plaquettes near a free boundary 
may be connected even when they are far apart from 
each other, increasing the number of possible connec- 
tions. On the other hand, full peb cannot be realized 
with the matching algorithm. Hence, we have limited 
our system sizes to L < 1024, M < 16 for R > 1, and 
L < 12, M < 768 for R < 1. 



III. RESULTS 

For Gaussian interactions we have considered systems 
with P-AP and F-DW boundary conditions with aspect 
ratios R — 1/8 to i? = 64. We took system sizes L = 1 to 
L = 12 for R = 1/8, 1/4, 1/2, and M = 1 to M = 16 for 
R = 1,2,4,..., 64. Typically ~ 10 5 samples were treated 
per system size (L, M), except for few combinations, no- 
tably R = 64 where we studied ~ 10 4 samples. 

For bimodal interactions we have studied P-AP bound- 
ary conditions for M = 1 to M = 16 (R = 1, 2, 4, . . . , 64), 
with ~ 10 4 samples per system size. 
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FIG. 1: Average domain-wall energy (Eow) of Gaussian sys- 
tem with P-AP boundary conditions as function of width M 
for different aspect ratios R. Lines are guides to the eyes only. 
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FIG. 2: Average domain-wall energy (Eow) of Gaussian sys- 
tem with F-DW boundary conditions as function of width M 
for different aspect ratios R. Lines are guides to the eyes only. 

We start with the results for the Gaussian distribution. 
In Figs. ^] and || the scaled average domain-wall energies 



L(E BW ) are displayed for P-AP and F-DW 



conditions respectively. From the scaling theorya a M 1+ 



boundary 
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behavior is expected for JJ > 1. Indeed for larger val- 
ues of R and M straight lines are visible in the double 
logarithmic plot. From fitting to algebraic functions, we 
have obtained for both cases the values of 9 as a func- 
tion of R; for small aspect-ratio the small system sizes 
were omitted from the fits. The results are displayed in 
Table |. 



R 


dp-AP 


#F-DW 


"F-DW 


1 


-0.285(3) 


-0.271(2) 


-0.153(2) 


2 


-0.288(3) 


-0.279(4) 


-0.215(2) 


4 


-0.288(1) 


-0.286(2) 


-0.249(2) 


8 


-0.284(2) 


-0.282(2) 


-0.265(2) 


16 


-0.288(2) 


-0.284(3) 


-0.273(2) 


32 


-0.289(2) 


-0.289(3) 


-0.274(3) 


64 


-0.290(7) 


-0.288(7) 





TABLE I: Resulting stiffness exponents from fits to functions 
L(Edw) = aM 1+e for P-AP and F-DW boundary conditions. 
The final column gives the result for F-DW boundary con- 
ditions when pbc are applied in the y-direction (data from 
Ref. ||). 

For the P-AP case, the exponent 9 is more or less in- 
dependent of the system size. By contrast, for the F-DW 
case a significant increase with system size is observed, 
with the value seeming to converge near the value of the 
P-AP case for laxge R. This i?-dependence explains why 
in previous workQ (R = 1) a difference between P-AP and 
F-DW has been found. The result found here is compat- 
ible with the exponents being equal for both cases. We 
quote 9 = —0.287(4), which is a conservative estimate 
including all values of 9 found for the P-AP case. Thk 
is compatible with the result of 9 = —0.282(3) obtainedEl 
for systems with full periodic be. 

Part of the data in Ref. ^| corresponds to the F-DW 
case with pbc in the y-direction, and we show this here 
in the final column in Table B where M < 12 was used. 
For this case an even stronger dependence of the effective 
exponent on R is observed, indicating that corrections to 
scaling are even larger in this case. We suggest that this 
is because the pbc forces the domain wall to take at the 
top and bottom edges the same positions along the x- 
axis, thereby raising its typical energy. In Ref. o it was 
shown that #p_ DW (called #f-af in Ref. ^ nevertheless 
extrapolates, for R — > oo, to an asymptotic value consis- 
tent with that obtained from the present work, where we 
employed free boundary conditions in the y-directian. 

We have also checked the scaling predictions^ for 
(£dw) explicitly. In Fig. || the scaled domain- wall en- 
ergy M~ 9 (Er,w) is plotted for all considered values of 
M as a function of the aspect ratio R. A very good data 
collapse can be observed (see e.g. right inset of Fig. || 
for P-AP). The functional form of the scaling function 
agrees very well with the predictions shown in Eq. ([}]). 

In Fig. [5| the normalized standard deviation of the scal- 
ing plot is shown as a function of R. Again, the high qual- 
ity of the scaling can be seen. Note that the statistics for 
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FIG. 3: Scaling plot of M- (E D w), with 9 = -0.287, as a 
function of R for P-AP (main plot) and F-DW (left inset) 
boundary conditions for all values of M which have been con- 
sidered. The curve is indeed a (very good) collapse of the 
results for different values of M, as exemplified for the P- 
AP case by the right inset, which shows a blow-up of the 
R = 2 points. The lines represent functions ~ R~ 1 (con- 
tinuous lines), ~ i? 9- ^ -1 )/ 2 (dashed line, main plot), and 
~ i? 9 -^- 1 ) (dashed line, left inset). 
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FIG. 4: Standard deviation a(R) of scaling plot points dev- 
ided by mean m(R). The mean was taken over all considered 
values of M. Lines are guides to the eyes only. 

both P-AP and F-DW cases are comparable, while the 
statistics for R = 64 are worse than for the other val- 
ues due to the smaller number of samples. From Fig. ^ 
we see that, for small values of R, the F-DW case scales 
less well than the P-AP case. This may explain, why, for 
F-DW, the stiffness exponent obtained for small aspect- 
ratios differs significantly from the R — > oo limit. The 
fact that F-DW scales less well than P-AP is probably 
due to the larger influence of the system boundaries in- 
duced by the free be. As a result, the scaling near R = 1 
is more complicated than a simple M 1+e behavior. 
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FIG. 5: Domain-wall energy (Bdw) of ± J system with P-AP 
boundary conditions as a function of L for M — 4. The line 
represents the function g(L) = 2.27 exp(— 0.59L). 

Finally, we turn to the system with bimodal distri- 
bution of the bonds. In Ref. || a saturation of the 
domain-wall energy has been observed for R — 1 when 
L = M —> oo, i.e. a fundamentally different result from 
the Gaussian case. This raises the question of whether 
the aspect-ratio approach is also able to make those re- 
sults consistent. 

As an example, in Fig. the domain-wall energy is 
plotted for M = 4 as a function of L. Clearly an exponen- 
tial decrease with L can be observed. This is due to the 
discrete nature of the interactions and can be explained 
as follows. For L » M, the system can be decomposed 
into almost independent subsystems of size M x M. Let p 
be the probability to find a non-zero domain- wall energy 
in a subsystem. Then the probability to find a non-zero 
J?dw in the full system is p L / M which decreases expo- 
nentially with L. 

As a result, for large R (which implies a large value of 
L for fixed M), the probability to find a non-zero domain- 
wall energy is so small that for a reasonable number of 
samples -Edw will be exactly zero. This is what we have 
indeed found for R = 32 and R = 64. It means that 
for the bimodal model the aspect-ratio approach cannot 
be applied. This is in contrast to the Gaussian case, 
where aspect ratio scaling works better and better with 



increasing R. 



IV. SUMMARY AND DISCUSSION 

We have studied the scaling behavior of the energy of 
domain walls in two-dimensional spin glasses obtained 
by changing the boundary conditions at T = 0. For 
the numerical calculations, we have employed a very effi- 
cient matching algorithm, enabling us to calculate exact 
ground states of systems up to size N = 1024 x 16. To 
eliminate finite-size effects, we have applied the aspect- 
ratio scaling approach. Our main results are: 

(i) The value of the stiffness exponent for Gaussian 
distribution of the bonds in d = 2 seems to be in- 
dependent of the choice of the boundary conditions 
and the way the domain walls are introduced. We 
give a final value of 9 = -0.287(4). 

(ii) The scaling behavior predicted by the aspect-ratio 
approach applies very well for the Gaussian distri- 
bution of the bonds, again indicating that the value 
of 9 is independent of the boundary conditions. 

(iii) The aspect-ratio technique does not allow us to gain 
further insight for systems with the bimodal distri- 
bution of the bonds. 
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